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ABSTRACT 

"<N, 

The empirical relationship between the broad line region size and the source 

^ ! luminosity in active galactic nuclei (AGNs) is used to obtain black holes (BH) 

masses for a large number of quasars in three samples. The largests BH masses 
iio 



Q , found exceed 10 M and are correlated, almost linearly, with the source lumi- 

nosity Such BH masses, when converted to galactic bulge mass and luminosity, 
indicate masses in excess of 10 13 M and a* in excess of 700 km/sec. Such massive 
galaxies have never been observed. The largest BHs reside, almost exclusively, 
in high redshift quasars. This, and the deduced BH masses, suggest that several 
scenarios of BH and galaxy formation are inconsistent with the observations. Ei- 
ther the observed size-L relationship in low luminosity AGNs does not extend 
to very high luminosity or else the M B h — M hu \ ge — a* correlations observed in 
the local universe do not reflect the relations of those quantities at the epoch of 
"^h galaxy formation. 

Q-f 

Subject headings: black hole physics — galaxies: active — galaxies: nuclei — 
to '. galaxies: high redshift — quasars: general 



1. Introduction 

Recent progress in reverberation mapping of active galacic nuclei (AGNs) allowed the 
first meaningful correlation between the broad line region (BLR) size (-Rblr) and the black 
hole (BH) mass in more than 30 objects. This provided a simple way to calculate BH masses 
for a large number of sources and resulted in a flood of papers on this topic. Some papers 
(e.g. Vestergaard 2002, hereafter V02; McLure and Jarvis 2002) investigated, in great detail, 
the wavelength dependence of the -Rblr — L — M relationship and provided useful ways for 
adopting the method to other wavelength bands. This opens the way for the study of BH 
masses in large samples of high luminosity high-z quasars. 
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All the new BH mass estimates are based on a single relationship obtained for a single 
sample of 34 AGNs for which BLR sizes are available from decade long reverberation mapping 
campaigns. More than half the sample was observed at the Wise observatory over a period 
of about 12 years (Kaspi et al. 2000, hereafter K00). Other objects have been monitored 
in other observatories and in several "AGN watch" campaigns (Netzer & Peterson 1997; 
Peterson 2001). The main findings are a significant -Rblr — ALa(5100) relationship (La(5100) 
is the monochromatic luminosity at 5 100 A) and the confirmation that the BLR gas is in virial 
motion (e.g. Peterson and Wandel 2000). These, plus the (model dependent) conversion of 
the observed full-width at half maximum (FWHM) of various emission lines into 3-D gas 
velocities, are sufficient to derive the mass of the central BH. 

This letter discusses the mass of the largest BH in the universe; those found in the 
centers of the most luminous quasars. It follows the works of Laor (1998; 2001), McLure 
& Jarvis (2002), Woo & Urry (2002) and others who used such methods for obtaining BH 
masses beyond the original K00 sample. The paper addresses also the Shields et al. (2002) 
new results and extends the mass estimates to much larger quasar samples. Section 2 presents 
new mass calculations for a large number of sources and §3 illustrates the new correlations 
found. Section 4 discusses the new results in light of the available information on the largest, 
most luminous elliptical galaxies and the epochs of quasars and galaxy formation. 



2. The largest BH 

2.1. BH Mass measurements 

New mass estimates have been obtained for a large number of AGNs using the -Rblr — L 
relationship obtained from the K00 sample; the only sample available for such calibration. 
This relationship is given, schematically, by 

#BLR = C X L\ (1) 

which results in the following mass estimate: 

M BH = c 2 L\[FWHM} 2 . (2) 

Here C\ and c 2 are constants that include the flux normalization and various assumptions 
about the velocity field in the BLR. The slope 7 is derived from the reverberation campaigns 
results and is in the range 0.5-0.7 (see below). The expression in eqn. 2 can be used to 
derived "single epoch" masses that combine the constants 7 and c 2 with observed FWHM of 
certain emission lines in individual objects. The method has been described in various papers 
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including KOO, V02, and McLure & Jarvis (2002). Its more useful applications are based on 
measured AL A (5100) and FWHM(H/3) for low redshift sources (the quantities used by KOO) 
and the combination of AL A (1350) and FWHM(C IV A1549) for high redshift objects. V02 
has looked into the inter-calibration of the two and supplied the expressions that are used in 
this work except for a small correction in the value of C2 that was introduced to adjust her 
constants to the cosmology assumed here: H Q = 70 kms -1 Mpc" 1 , Q m = 0.3 and fl\ = 0.7. 
McLure & Jarvis (2002) provided similar expressions for Mgll A2798 which are not used in 
this work. 



2.2. The sample 

Three AGN samples have been used in this work: 1) the LBQS sample (Forster et al. 
2000 and references therein), 2) A sample of 104 high redshift high luminosity quasars with 
ground-based spectrophotometry (L A (1350)) and good FWHM(C IV A1549) measurements, 
and 3) the new L A (5100) and FWHM(H/5) listed by Shields et al (2002). Many of the 
sources in the second sample are UM quasars and the raw data can be found in MacAlpine 
and Feldman 1982, Baldwin, Wampler & Gaskell (1987) and Baldwin (1977). 

Forster et al. (2000) supplied monochromatic luminosities and FWHMs for many emis- 
sion lines in about 1000 LBQS quasars. Since many of the sources have been observed 
through relatively small aperture, and under poor weather conditions, it was decided to use 
the Bj magnitudes that are much more accurate. This follows Green et al. (2001) who 
studied the Baldwin relationship in this sample and obtained monochromatic luminosities 
using the same method. All fluxes have been corrected for galactic reddening using the 
Green et al. procedure. A major assumption here, and in Green et al. (2001), is that the 
observed continuum can be described by a single L u oc v~ a power-law with a = 0.5. This 
approximation neglects the possible dependence of a on source luminosity which may affect 
the L — M relationship (see §4). Forster et al. (2001) provided several different measure- 
ments of FWHM(C IV A1549) with and without the narrow line component. The "single" 
component fit was used and the "broad only" fit was checked to verify that the results 
are not sensitive to this choice. A handful of sources with FWHM(H/3)< 1,000 km s _1 or 
with FWHM(C IV A1549)> 20,000 km s _1 were removed from the sample since those were 
considered unreliable or affected too much by the narrow emission line. As for the second 
C IV A1549 sample, no galactic reddening was applied and the same assumption about L v 
was used. In this case there is no significant dependence on a since the original papers quote 
the observed flux at around rest wavelength of 1450A. 

The above samples are optically selected and suffer from various selection effects. This 
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is of no real consequence to the main goal of the paper which is to derive the mass of the 
largest known BHs. It may affect, however, the derived M — L correlations (§4). 

3. The L — M relationship for high luminosity AGNs 

BH masses have been calculated using equation 2 and the normalizations derived by 
K00 and V02 adjusted to the cosmology chosen here. The determination of the slope 7 is 
crucial for the present work and will be discussed prior to presentation of the new results. 

We start from the original K00 sample to which we apply two statistical methods for 
finding 7: the Akritas & Bershady (1996) BCES estimator (for which we only consider 
the BCES bisector) and the fitexy method described in Press et al. (1992). The merits 
of the different methods have been discussed, extensively, in several papers and will not be 
repeated here. Our experience shows that that the differences between the slopes obtained by 
the different methods are larger than the formal uncertainties on the slopes of each method. 
The K00 sample adjusted to the new cosmology gives 7 = 0.58 ±0.12 for the BCES bisector 
estimator and 7 = 0.68 ± 0.03 for the fitexy method. The two are formally consistent with 
each other and 7(BCES) was adopted here. 

Since the purpose of this work is to extrapolate to very large L, we also experimented 
with removing the lowest luminosity objects from the sample. Removing the three objects 
with AL A (5100) < 10 43 ergs s" 1 resulted in 7 (BCES) = 0.71 ± 0.21 and f(fitexy) = 0.69 ± 
0.03. Removing the seven objects with AL A (5100) < 10 43 7 ergs s" 1 resulted in 7(BCES) = 
0.58 ± 0.19 and ^(fitexy) = 0.74 ± 0.04. All these results suggest that the two methods are 
consistent with each other and the slope cannot be determined to an accuracy better than 
about 0.15. The value adopted for illustrating the results of this work is the smaller one 
found for the entire K00, 7 = 0.58. The implications for the case of larger or smaller 7 are 
discussed in §4. 

Shields et al. (2002) suggested the use of the "physically motivated" value of 7 = 0.5. 
The strongest argument for using this value is the suggestion by Netzer & Laor (1993) 
that the outer boundary of the BLR is determined by the dust sublimation radius which 
is similar to the measure -Rblr to within a factor ~ 2. There are several problems in 
applying this idea to the present mass determination. First, the "reverberation radius" is 
determined by the responsivity of H/3 to changes in the ionizing luminosity, L ion , which is 
smaller than the bolometric luminosity that determines the dust sublimation radius. In 
addition, 7 = 0.5 means the same BLR ionization parameter for low luminosity Seyferts 
and the highest luminosity quasars. This has never been shown to be the case in large QSO 
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samples. Thus, more work is required to justify this theoretical value of 7. 

The masses computed with the 7 = 0.58 slope are presented in Fig. 1. The diagram 
contains data for 505 QSOs with CivA1549 measurements and 219 source with H/3 mea- 
suremets. The luminosity range is roughly AL A (1350) = io 44 ~ 47 - 5 ergs s -1 . Also shown is 
the best regression line (see below) and the mass range of ±a M around the median calcu- 
lated in luminosity bins of 0.3 dex. The largest BHs are found in sources with z > 2 with 
M B h — 10 10 ' 2 M (15 with mass exceeding 10 10 M Q ). Using 7 = 0.68 (the slope found with 
the fitexy method), raise this number to about 10 104 M (62 with mass exceeding 10 10 
M Q ). 




XLJ1350 A) 



Fig. 1. — Black hole mass as a function of ALa(1350) for the quasar samples described in the 
text assuming 7 = 0.58. Open symbols represent mass obtained from the H/5 line and full 
symbols masses obtained using CivA1549. The two CivA1549 samples completely overlap 
in properties and were not given different symbols. The dashed lines represent the ±cr M 
range around the median and the straight line the M cx L 8 relationship. 
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Fig. 2 shows Mbh as a function of redshift for the same sample under the same assump- 
tions. 




z 



Fig. 2. — Mbh vs. redshift for the sample in figure 1 with the same symbols. 



The data in Fig. 1 suggest a simple linear dependence of the form Mbh oc LP . This 
has been tested by performing a linear regression analysis using the same two methods 
described earlier. The procedure used for calculating the errors is the following: For L1350, the 
assumption is of a constant error of 0.15 dex representing the measurement uncertainty, the 
extrapolation in wavelength and the typical range in luminosity due to continuum variability. 
This number does not affect the resulting slope (5 in any significant way. As for the mass, this 
was done using standard error propagation combining all errors due to the uncertainties in L 
and in FWHM (line width uncertainties are given in Forster et al. 2000). The combined error 
for this case is typically 0.15-0.25 dex. No uncertainties are listed for FWHM(C IV A1549) in 
the second quasar sample and for FWHM(H/3) in Shields et al. (2002). A uniform error of 
0.2 dex in M B h was assumed in those cases. The errors are relatively large and are expressed 
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in logarithmic form (i.e. 0.5(log(x+dx)— log(x— dx), see Lyons, 1991). 

Table 1 lists several slopes obtained by the two methods for our standard case of 7 = 0.58 
and for 7 = 0.68, the slope obtained by the fitexy method. Given the various biases and 
unknowns, it is reasonable to assume that the real uncertainty in (5 is at least as large as 
±0.15. With this uncertainty, the slopes of the CivA1549 sample and the entire sample are 
barely consistent with each other and the slopes of the H/5 sample and the entire sample 
are indistinguishable. The scatter in slope is probably due to the very different luminosity 
range of the C IV A1549 and the H(3 samples. A second approach that was tries assumed a 
uniform uncertainty in Mbh of 0.3 dex for all objects. This gave very similar results. The 
overall conclusion is that for luminous AGNs, Mbh oc L 9±015 . 

The Mbh — L correlation found here is very different from the one found in K00. The 
reason is probably the incompleteness of the small K00 sample which resulted in a biased 
sampling of FWHM(H,3) vs. ALa(5100) not representing the parent population. Indeed, K00 
found FWHM(H/9)oc L~°- 27 while in the samples under study the correlation is much flatter. 
The FWHM-luminosity dependence in various samples will be addressed in a separate paper 
(Corbett et al. 2003). 

The tight Mbh — L relationship enables the study of the Eddington ratio, L/L^dd, in 
these samples. The observed M — L relationship suggests a very weak, if any dependence 
of L/L Edd on L or on M B h- This impression is confirmed by a formal statistical analysis. 
Since the results are marginal, they will not be presented here. Another important issue 
is the mean L/L-^dd- This depends on the distribution in this property as well as on the 
exact conversion from \L\ to bolometric luminosity and the value of 7. Assuming first 
7 = 0.58 and L^dd = 9ALa(5100), as in K00, gives a median L/L^dd of 0.53. The composite 
spectra published recently by Telfer et al. (2002) suggest a different conversion with L^d — 
5ALa(5100). This translates to a median of 0.28. The above values are transformed to 0.33 
and 0.18 for the case of 7 = 0.68. In both cases the distribution wis wide, covering about 
a factor 10 in L^dd- Thus, the choice of 7 = 0.58 results in a large mean L/L^dd and a 
large number of sources with super Eddington luminosities. As explained by Woo and Urry 
(2002), the implications to the derived M B h — L relation are very importance (see §4). 



4. Discussion: the largest BHs and the most luminous galaxies 

The new results presented here suggest that the largest BHs are situated in the most 
luminous quasars that are, typically, the highest redshift sources. At the extreme end of the 
distribution we find BH masses of order 5 x 10 10 M if 7 = 0.7, and 1.1 x 10 10 M if 7 = 0.5. 
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This is greater than obtained so far in large samples. The recent work by Shields et al. 
(2002) aimed at the calibration of the the [O ill] A5007 line width as a bulge mass estimator. 
The method is based on the close agreement between FWHM([0 ill] A5007) and the stellar 
velocity dispersion a* at low luminosity and the K00 mass estimates at higher luminosity. 
Using this method and 7 = 0.5 (their Table 2) they find one object with Mbh exceeding 10 10 
M Q and several others approaching this mass. As shown in Fig. 1, the CivA1549 samples 
includes many more sources with such large masses. 

Before addressing the cosmological consequences we note the various factors influencing 
the M — L relationship and likely reasons for overestimating Mbh- 

1. The K00 sample covers a limited luminosity range and all mass estimates corresponding 
to ALa(1350) > 10 46 are necessarily obtained by extrapolation. Since this is the only 
sample available so far, there is no independent way to verify the largest masses until 
successful reverberation mappings are obtained for higher luminosity AGN. Moreover, 
as explained in §2, the slope of the -Rblr--^ relationship is uncertain. The slope chosen 
here (7 = 0.58) is close to the middle of the range. Its increase to 0.7 will increase the 
mass at the high luminosity end by a factor of about 2.5. 

2. The largest new mass estimates are based on the measured AL A (1350) which is scaled 
to the K00 luminosity assuming the same spectral energy distribution (SED) for high 
and low luminosity AGNs. This assumption has never been tested in large quasar 
samples. The data for such test are already available (Telfer et al. 2002) but the 
results are not yet known. Intrinsic reddening, in the quasar host galaxy, is another 
potential complication related to the inter-calibration of optical and UV luminosities. 

3. The FWHM(CivA1549) may not reflect the virial motion of the BLR gas in high 
luminosity quasars. 

4. The samples used here suffer from various selection effects. This influence only slightly 
the largest derived masses but can affect much more the L — M B h correlation. For 
example, magnitude limited samples may not include the less luminous quasars, those 
with the smallest L/L Edd . This results in a false impression of a very strong M — L 
correlation. Woo and Urry (2002) investigated this idea in great detail and concluded 
that all strong M — L correlations obtained so far suffer from such a selection effect. 

The main conclusion of this work is that the largest BH masses are found in the highest 
luminosity quasars. The masses of such BHs can reach the extreme values of io 10 - 3 ~ 10 - 6 M , 
depending on the value of 7. Using recent conversions to host galaxy properties one finds 
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M bu i ge ~ io 13 - 1 - 13 - 4 M (Kormendy & Gebhardt 2001), M B ,bui g c ~ -25 mag (Kormendy 
2001) and a* exceeding 800 km/sec (Tremaine et al. 2002). Such galaxies have never been 
observed and are not predicted to exist by standard galaxy formation theories. 

In principle, this is still consistent with the observations since the sources with the largest 
Mbh are the most luminous ones and will completely out-shine any host galaxy. Thus, there is 
no direct way to rule out the existence of such galaxies. However, the theoretical implication 
are in conflict with recent ideas that the largest galaxies attain their mass through a series 
of mergers, a process that operates continuously to redshift 2 or smaller. A similar difficulty 
is found for the BH growth since those same theories (e.g. Haehnelt and Kauffman 2000; 
Yu and Tremaine 2002) assume that galactic nuclei BHs increase their mass up to redshifts 
smaller than 2 by the same series of mergers (or, perhaps, only through large mergers). Thus, 
the largest BHs are predicted to be associated with the most massive galaxies at z < 2, in 
conflict with the data in Fig. 2. It is clear that active BH with M B h > 10 10 M are not 
found in the local universe. It is also clear that dormant BHs of this mass, or the galaxies 
with extreme properties that are supposed to host such BHs, have never been found. The 
whereabout of the huge BH formed at z ~ 3 is thus unknown. 

A more plausible suggestion is that some or all the conversion factors used to obtain 
the galactic mass, magnitude and a* from the BH mass, that are based on measurements in 
the local universe, cannot be extrapolated to high luminosity high redshift objects. Perhaps 
they are only valid at z < 2, after galaxies and nuclear BHs have accumulated most of their 
mass. If correct, this would mean that some "normal looking" galaxies contain extremely 
massive BHs. A similar suggestion by Laor (2001) involves a dependence of Mbh/ ^tmige on 
the BH mass or the absolute magnitude of the host galaxy. 

To conclude, either the measurements of BH masses presented here for the most lumi- 
nous quasars are grossly overestimated, because of the reasons described above, or else the 
relationships between BH masses and various properties of their host galaxies at high z are 
very different from those measured in the local universe. A second conclusion, which is less 
certain because of various selection effects, is that for AGNs, M BH oc £°- 9±0 - 15 . 

The work described in this paper is based primarily on a decade-long AGN monitoring 
project at the Wise observatory. I am grateful to many of my colleagues and students 
that helped in making this into a very successful project. Special thanks go to Dan Maoz, 
Shai Kaspi and Ohad Shemmer who led various parts of the project and without whom it 
would have been impossible to bring it to completion. Useful discussions with Ari Laor are 
gratefully acknowledged. This work is supported by the Israel Science Foundation grant 
545/00. 
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Table 1. Regression analysis results for Mbh = aL 



method 


7 


sample 


P 


a 


BCES Bisector 


0.58 


all C iv A1549 


1.12±0.05 


-42.9 


fitexy 


0.58 


all Civ A 1549 


1.13±0.04 


-43.1 


BCES Bisector 


0.58 


all H/3 


0.83±0.03 


-29.2 


fitexy 


0.58 


all R(3 


0.80±0.02 


-27.6 


BCES Bisector 


0.58 


all objects 


0.80±0.02 


-27.7 


fitexy 


0.58 


all objects 


0.78±0.01 


-26.9 


BCES Bisector 


0.68 


all objects 


0.90±0.02 


-32.2 


fitexy 


0.68 


all objects 


0.89±0.01 


-31.9 



